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We present a microscopic theory of the magnetic field induced mixing of heavy-hole states ±3/2 in 
GaAs droplet dots grown on (111) A substrates. The proposed theoretical model takes into account 
the striking dot shape with trigonal symmetry revealed in atomic force microscopy. Our calcula- 
tions of the hole states are carried out within the Luttinger Hamiltonian formalism, supplemented 
with allowance for the triangularity of the confining potential. They are in quantitative agreement 
with the experimentally observed polarization selection rules, emission line intensities and energy 
splittings in both longitudinal and transverse magnetic fields for neutral and charged excitons in all 
measured single dots. 

PACS numbers: 72.25.Fe, 73.21. La, 78.55.Cr, 78.67.Hc 



I. INTRODUCTION 

Thorough understanding and tuning of the optical 
properties of semiconductor quantum dots (QDs) has al- 
lowed the development of efficient quantum emitters)^ In 
parallel, the electrical and optical manipulation of the 
spin states of single carriers in QDs has made unprece- 
dented control ove r quantum states in the condensed 
matter possible Single spin manipulations and quan- 
tum optics with single photons are directly related to the 
optical selection rules. 4 The band structure given by the 
symmetry and strength of the carrier confinement po- 
tentials defined, e.g., by the dot shape and its chemical 
composition plays a crucial role in defining the exact car- 
rier spin state and, correspondingly, the polarization of 
the emitted and absorbed photons. 




FIG. 1: (Color online): AFM images of GaAs quantum dots 
grown on (111) A substrates reveal the triagonal symmetry of 
the dots, with shapes varying from irregular hexagons to equi- 
lateral triangles depending on crystallisation temperature. 

The majority of devices based on self-assembled quan- 
tum dots is grown along the standard [001] crystal axis, 



for which the dots tend to have an elliptic or rectangular 
shaped base. The associated C2v point symmetry results 
in a complex exciton fine structure^^ characterized by 
well-defined linearly polarized eigenstates being split in 
the case of neutral X° exciton emission by the electron- 
hole exchange interaction. The main aim here is the sup- 
pression of this fine structure splitting to allow the emis- 
sion of entangled photon pairs, as convincingly demon- 
strated using a variety of techniques.®^ Such a low, C2„ 
point symmetry is also relevant for the [001]-grown quan- 
tum dots with circular or square-shaped base provided 
that the top and bottom interfac es ar e non-equivalent 
or if an in-plane strain is present .12112103] A promising, 
alternative approach to suppress the fine structure split- 
ting of the neutral exciton is the growth along the [111] 
crystal axis, which is also the orientation of most nano- 
wiresP^This growth axis has the advantage of providing 
in principle quantum-dot shape of the higher point 
symmetry. Really, small fine structure splittings in as 
grown [1 11] q uantum dot structures have been recently 
predicte d 16 * 17 ^ and observecP^HH] making such structures 
a very promising system for entangled photon pair emis- 
sion. 

The exact electron spin states, selection rules and op- 
tical properties in any semiconductor nanostructure are 
primarily determined by the symmetry of the investi- 
gated systemPSl As compared to the well studied [001] 
grown quantum dots, the triagonal symmetry associated 
to the growth along the [111] axis revealed in the Atomic 
Force Microscopy (AFM) images presented in Fig. [T] 
gives rise to surprising new spin effects discovered re- 
cently and studied in detail in the present paper: 

(i) In a magnetic field applied along the [111] growth 
axis (Faraday geometry), the neutral and singly charged 
exciton complexes show a total of four interband transi- 
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tions each, two a + and two a~ polarized,^ see Fig. [2J1 
and b, left panels. 

(ii) In a transverse magnetic field applied perpendic- 
ularly to the growth axis (Voigt geometry), the charged 
exciton emission is split into two lines, each doubly de- 
generate, see Fig. [2^ and b, right panels. Experiment 
shows that the heavy hole states have zero Zeeman split- 
ting in the transverse field. 

Both of these findings, (i) and (ii), detailed in section [TT| 
are in stark contrast to results obtained in samples grown 
along the [001] axis (see, for example, Refs. I5|25l - l27|l . 
where two transition lines are observed in the Faraday, 
and four well separated transitions in the Voigt geometry. 



Symmetry arguments, presented in Sec. Ill A allow us to 
identify the key ingredients needed to explain the char- 
acteristic interband transitions: the valence heavy hole 
pseudo spins | +3/2) and | — 3/2) are coupled in longi- 
tudinal magnetic field, whereas they remain uncoupled 
in transverse field. Hence, the magnetic field induced 
tuning of the energy level structure and the underlying 
symmetry in [111] grown quantum dots, see Fig. [2J; and 
d, makes it possible to develop new, original approaches 
to coherent electron and hole spin controlJ^smfll 

The main part of this paper, Sec. |III B| is devoted to 
a detailed microscopic model that provides quantitative 
evaluation of the magnetic field induced heavy-hole mix- 
ing in the Faraday geometry. We show that the main 
contribution to the heavy-hole mixing results from the 
influence of the triangular quantum dot shape on the en- 
velope part of the valence states wave- functions, rather 
than the underlying lattice symmetry due to the [111] 
growth axis. Details of the heavy hole Zeeman effect cal- 
culations are provided in the Appendix. 

The experimental and theoretical results presented so 
far open up for the true quantum mechanical engineering 
of the Zeeman effects in nanostructures. 



II. EXPERIMENTAL RESULTS 

The samples were grown by droplet epita xy usin g a 
conventional molecular beam epitaxy syste m 18 ! 31 ! 32 ! on 
a GaAs(lll)A substrate. The GaAs dots are grown on 
100 nm thick Alo.3Gao.7As barriers and are covered by 
50 nm of the same material. The barrier and the quan- 
tum dot material have the same lattice parameter} 3 ^ the 
samples are within good approximation strain free. As 
a result, there are no piezo-electric effects to take into 
account, despite the [111] growth direction. There is 
no continuous wetting layer in the sample connecting 
the dots, whose typical height is ~ 3 nm, and radius 
is ~ 15 nm, see details in Ref. I18[ which results in strong 
confinement of the carriers in the dots. Atomic force 
microscopy studies carried out on similar samples grown 
on a GaAs(lll)A substrate reveal a dot base with trig- 
onal symmetry, see Fig. [I] that evolves from irregular 
hexagons to equilateral triangles as the crystallization 
temperature is varied.^ Single dot photoluminescence 
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FIG. 2: (Color online) (a) QD1 Left: Measured transition 
energies as a function of the applied magnetic field B || [111]. 
Right: Measurements for B _L [111], (b) Same as (a), but for 
QD5. (c) Polarization selection rules for X + recombination in 
the Faraday geometry (B || [111]) with mixed hole states, see 
Eq. Q, where Ci,2 are coefficients determining hole eigen- 
states. (d) Polarization selection rules for X + recombination 
in the Voigt geometry (B _L [111]). 



(PL) at 4 K is recorded with a home build confocal 
microscope with a detection spot diameter of ~ 1 (xm. 
The detected PL polarization is analysed, the signal is 
dispersed by a spectrometer and detected by a Si-CCD 
camera. Optical excitation is achieved by pumping the 
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AlGaAs barrier with a HeNe laser at 1.96 eV that is lin- 
early polarized to exclude the effects of optical carrier 
orientation and dynamic nuclear polarization! 35 ! 36 ! 

Three exciton complexes dominate the PL emission at 
zero applied field: The already mentioned neutral X°, 
the negatively charged exciton (trion) X~ comprised of 2 
electrons in a spin singlet state and 1 hole, and the posi- 
tively charged trion X + (1 electron, 2 holes in a spin sin- 
glet state) are identified using fine structure analysis and 
optical orientation experiments The high symmetry of 
the dots grown along [111] is reflected in typical values 
for the splitting of the X° emission due to anisotropic 
electron-hole exchange Si of a few //eVj^ extracted from 
angle dependent PL polarization analysis in the linear 
basis. 

The X + and X~ trions of all dots investigated exhibit 
four transitions in longitudinal magnetic fields, as seen 
from the emission spectra of two typical dots presented 
in Fig. [2] Surprisingly, there are clear polarization se- 
lection rules: two transitions are o~ polarized, the two 
others a + . The competition between the exchange inter- 
action and magnetic field effects leads to intriguing obser- 
vations for the neutral exciton X°: The Zeeman splitting 
is strongly non- monotonous,^ see Fig. [3] for a striking 
example. 

Figure [2] shows also the spectra of the same two typical 
QDs in the Voigt geometry: For QD1 shown in panel (a) 
the total number of lines we are able to resolve is two 
for each charged exciton (X + and X - ), and not four as 
could be expected by analogy with [001] samples! 5 ! 25 " 27 ! 
This is true both in the circular (<j + /a~) and the linear 
basis (ir x /ir v ). Analysing a quantum dot with lower sym- 
metry allowed us to uncover the true polarization basis 
and lift the remaining degeneracies. Indeed, the shape 
of QD5 probably deviates from an ideal pyramid with 
an equilateral triangular base and each line of the ini- 
tially observed doublets undergoes a further very small 
splitting with the increasing field. According to Fig. [2]d 
at 9 T we observe four resolution limited transitions in 
the linear basis for the X + and X - trions with the ad- 
ditional splitting AE max ss 40 fieV. As shown below the 
presence of four transitions in the Voigt geometry indeed 
corresponds to the non-zero in-plane heavy hole g-factor 
allowed for (111) structures with point symmetry lower 
than C$ v . 

Moreover, both in the Faraday and the Voigt geome- 
tries nominally dark exciton (with total spin projection 
m z = s z + j z = ±2, where s z and j z are the electron 
and hole spin components) gains an oscillator strength, 
allowing a precise determination of the dark-bright X° 
splitting in the range of 300 to 400 fieV (determined for 
several tens of QDs), which is strongly enhanced, due to 
confinement, compared to the GaAs bulk value and the 
larger GaAs interface fluctuations dotsP^I 



III. THEORY 

A. Symmetry Analysis 

From the analysis of the single dot PL in linear basis, 
see Ref. [35] for details, we are able to extract an X° fine 
structure (anisotropic) splitting Si of typically only a few 
/ieV. This, together with the analysis of possible symme- 
tries for zinc-blende lattice based structures grown along 
z || [111] axis, allows us to conclude that the investigated 
QDs have (close to ideal) trigonal C% v point symmetry, 
in perfect agreement with the AFM measurements shown 
in Fig. ^ 

In what follows we introduce the Cartesian coordinate 
frame x || [112], y || [110] and z || [111] relevant for 
the system under study and analyze the effect of a lon- 
gitudinal magnetic field along the growth axis z. For a 
nanostructure of the C3„ symmetry we assume one of the 
planes a v £ C3 V to contain the axes x and z. The phases 
of heavy-hole basis functions |3/2), | — 3/2) are chosen to 
show the same point-group symmetry properties as the 
functions — j" (x + iy) and J, (x — iy) , respectively, where 
t, J, are the spin-up and spin-down columns. 

For dots with the C^ v symmetry, the group theory 
predicts a magnetic field induced mixing between the 
heavy hole states |3/2) and | — 3/2). These two states 
form a reducible representation reduced to the two one- 
dimensional representations T 5 and T 6 with the basis 
states |r 5 ) = |3/2}+i|-3/2) and |r 6 ) = |3/2) -i| -3/2). 
The latter functions are related by the time inversion op- 
eration and, therefore, at zero magnetic field have the 
same energy. Since the direct product 



(r 5 + r 6 ) x (r* + q) - 2Ti + 2r 2 



(i) 



contains two representations T2 according to which trans- 
forms the magnetic field component B Zl the heavy-hole 
Zeeman splitting in the basis |3/2), | — 3/2) is described 
by a 2x2 matrix with two linearly independent parame- 
ters 



Hi 



1 



VbB z 



9 hi 

9h2 



9h2 
-9 hi 



(2) 



Here \xb is the Bohr magneton, g^i and gya are the effec- 
tive hole g- factors, the latter also being real if the trigo- 
nal symmetry includes the mirror reflection in the plane 
(x,z) I] (110). The similar result holds for the trigo- 
nal point groups D3 and D^d = D3 x C,. For the low- 
symmetry trigonal C3 and Sq = C3 x C$ point groups, 
the mixing of the heavy holes by the longitudinal field 
is allowed as well. In this case, however, the symmetry 
does not fix the orientation of lateral axes x and y where 
the off-diagonal matrix elements are real. For all hexago- 
nal groups including C^h and D 3 ^ the longitudinal field- 
induced mixing of the heavy hole states is forbidden, and 
the value of gh2 in Eq. Q vanishes. This is also true for 
the orthorhombic C2V group and tetragonal T)2d group 
which reproduce the symmetry of the conventional [001] 
grown dots (lens- and disk-shaped, respectively). 
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TABLE I: (?-factors (typ. error < 10%) for charged and 
neutral excitons obtained from fitting the data. For the X° 
the g e and gh values obtained for X + for the same dot are 
taken and only \ghi\ is varied to fit the bright and dark X° 
splitting simultaneously. 





QD1 




QD5 




geometry 


Faraday 


Voigt 


Faraday 


Voigt 


X~ : g e 


0.48 


0.80 


0.49 


0.73 


9h 


0.79 





0.83 


0.07 


\9h2\ 


0.57 




0.53 




X+ : g e 


0.47 


0.84 


0.47 


0.79 


9h 


0.72 


0.03 


0.71 


0.07 


\9h2\ 


0.7 




0.62 




X° : |ffe| 


0.47 


0.80 


0.47 


0.77 


9h 


0.72 





0.71 


0.07 


\9h2\ 


0.59 




0.5 





For completeness we note that the experiments in the 
Voigt geometry shown in the right panel of Fig. [2^i and 
b infer that heavy holes are not mixed by the in-plane 
magnetic field. Indeed, the direct product (JIJ does not 
contain the representation T 3 describing the symmetry 
properties of the field components B x and B y and, there- 
fore, the point group does not allow for the mixing 
of the ±3/2 states by the in-plane magnetic field. This 
result is as well valid for all 12 point groups of the hexag- 
onal crystal family including the hexagonal and trigonal 
crystal systems. On the other hand, the mixing is al- 
lowed by the C2 V or D 2 d point-group symmetry where 
small but reliably measured in-plane heavy-hole g-factors 
were observed!^ 



1. Charged excitons 

In a longitudinal magnetic field, the hole eigenener- 
gies are E± = ±g h fi B B z /2 with g h = y/gf^ + g\ 2 and 
the hole eigenstates \h, ±) are admixtures of |3/2) and 
I — 3/2), as indicated in Fig.[2j:, with the coefficients 
determined solely by the ratio ghijghx, see Ref. l24l for 
details. For non-zero g^, all the four radiative transi- 
tions are allowed, each transition being circularly polar- 
ized, either a + or a~ . The transition energies include a 
Zeeman contribution as well as a diamagnetic shift com- 
mon for all lines. The splitting between different spin 
states is determined by combinations of the electron and 
hole effective g- factors whi ch allows to find a pair of pa- 
rameters, g e and gh = \/ 9hi + 9h2- The intensities of 
circularly-polarized lines are proportional to C\ and Cf 
and independent of the magnetic field, in full agreement 
with our experiments. From the ratio of intensities of 
identically-polarized lines we can find the ratio gm/gh 
and, therefore, determine values of g ei g^i and the abso- 
lute value of gh2^ 



2. Neutral exciton X° 

When analysing the emission energies and polariza- 
tion of the neutral exciton X°, a competition between 
the magnetic field and the electron-hole exchange inter- 
action has to be considered. Taking into account that 
in the absence of magnetic field the bright exciton states 
with the angular momentum components m z = ±1 are 
separated from the dark states with \m z \ = 2 by th e en- 
ergy So, we obtain for the X° sublevel energies^l2il 

E s z J = s z9c^bB z + -(6 +jS s J , (3) 

S Sz = \J$q + {ghHBB z ) 2 - 4:S z g h ifi B B z S . 

Here j = ±1 labels heavy-hole eigenstates in the 
presence of magnetic field and <5o > 0. It follows 
from Eq. (l3| that the splitting of bright X° states, 
-E+i/2,+i(-E?7) — E_i/2.+\{B z ) 1 can be a non-monotonous 
and even sign-changing function of B z as has been con- 
firmed experimentally!^ This result is another striking 
difference when comparing with the [001] grown dots, 
where the observed splitting increases mo noto nously as 
a function of the applied longitudinal field.^H 

As a particular example we present in Fig. [3] the tran- 
sition energies of QD6 determined in the Faraday geome- 
try. This dot shows comparatively weak heavy hole mix- 
ing, resulting in a relatively small value of gh2 — 0.42 
compared to gh = 1.16, and allows us to illustrate how 
well Eq. ([3]) fits the experiments: We even observe an 
anti-crossing of dark and bright X° states, see vertical 
dashed line in Fig. |3ji. Here the consequences of intro- 
duction of gh2 into the Zeeman Hamiltoninan, Eq. Q, 
and, hence, into the fit are striking. First, if gh2 is set 
artificially to zero, there is no anti-crossing but a clear 
crossing. This crossing appears at the characteristic field 
B z — So/(ghiHB)- The dark-bright separation at this 
field B z is exactly | gh2l^BB z \. Figure shows directly 
the evolution of the dark and bright transitions in 
the presence of heavy hole coupling, i.e., of non-zero gh2- 
Second, the small value of <?/ l2 leads to a strongly non- 
monotonous bright exciton splitting going from down 
to -150 /xeV at 6.5 T back up to -100 fieV at 9 T. The 
applied magnetic fields are too low h ere to see a reversal 
of the sign, but the tendency is clear^21i2l ( we observed, 
however, this sign reversal in other dots, see Ref. |2~4"}| . 
Hence, the high quality of the fits of the neutral exciton 
splitting in applied fields as presented in Fig. [3] under lines 
the important role of heavy hole coupling, i.e., non-zero 

.9^2 ■ 

It is worth stressing that all sets of experimental data 
obtained in the Faraday geometry for several tens of 
quantum dots, both for trion and exciton emission, is 
well reproduced within a simple model, where linear- 
in-B z contributions to the Zeeman effect with ghi and 
9h2 are taken into account onlyP^^l Thus, an in- 
troduction of only one additional parameter gh2 allows 
us to fit all the transition energy splittings of the neu- 
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subsequent estimation has shown that the correspond- 
ing constant *4( 3 ) is too small to be responsible for 
the effect. Here we propose a mechanism of the heavy 
hole mixing which is based on the standard Luttinger 
Hamiltonian quadratic in k and an anisotropic shape of 
the GaAs/AlGaAs(lll) QDs including polar orientation 
along the growth direction z and trigonal symmetry in 
the plane perpendicular to z. 



1. Quantization of the hole motion in the pyramidal 
quantum dot 

The hole states in a QD are described in the framework 
of the Luttinger Hamiltonian Hr 8 + V(r), where V(r) is 
the operator of potential energy of the hole and rir s is 
written in the standard matrix forirP 



Hi 



( F 

H* 
I* 

V o 



H 

G 


I* 



I 


G 

-H* 



\ 
I 

-H 
F J 



(4) 



In the spherical approximation used hereinafter, the func- 
tions F, G, H, and / read 



F={A~B)k 2 z 



A 



\k . 



FIG. 3: (Color online): QD6 (a) Transition energy of the 
nominally dark (squares) and bright (circles) neutral exciton 
X° state as a function of the applied magnetic field (Faraday 
geometry). Left- and right-hand sides represent the states 
with s z = 1/2 and s 2 = —1/2, respectively. Fit (green curve) 
of the data using Eq. [3] allows to extract So — 422 fieV, 
g h = 1.16, g e = 0.49, g hl = 1.05 and g h2 = 0.42. The fit in- 
cluded the diamagnetic shift for both curves as E = Eq + olB% 
with E = 1.819676 eV and a = 4.15 fieV /T 2 . Blue curve: 
identical parameters apart from gi^ = 0. (b) X° exciton 
bright splitting as a function of the magnetic field. Solid line 
is the fit obtained by using Eq. [3] 

tral and charged excitons, i.e., explain all the differences 
between the excitonic transitions observed in the [111] 
grown droplet dots and standard dots grown along [001]. 
In the symmetric [111] grown dots under study we of- 
ten find strong heavy hole mixing as \gh2\ > \ghi\- Such 
surprisingly high values of the off-diagonal g-factor are 
consistent with the microscopic theory developed in the 
next Subsection. 



B. Microscopical theory 

In Ref. [^| in order to explain a remarkable value of 
gh2 observed experimentally we surmised that the ±3/2 
heavy hole mixing could occur due to a cubic-in-fc con- 
tribution T-L V 3 to the hole Hamiltonian. However, the 



G=(A + B)k 2 z +(A-^\ (kl + kl), 

H = -V3Bk z (k x - iky) , / = — —B (k x — \k y ) . 

Here kj (j = x, y, z) are the Cartesian components 
of wavevector (we remind that the axes frame used is 
x || [112], y || [1 10] and z || [111]). We will further use 
hole representation so that constants A and B are re- 
lated to the Luttinger parameters 71 and 72 = 73 = 7 as 
A = ^ 2 7i/2to , B = ff^/mo, m o being the free electron 
mass. In the effective Hamiltonian Q the wavevector k 
is replaced by the operator — iV. 

We take the potential in the following matrix form 

(V hh {r) \ 

v()= V lh (r) 

V lh (r) 

V V hh (r)j 

allowing for different potential energies of the heavy (Vhh) 
and light (V^) holes. In self-organized QDs the size quan- 
tization along the growth axis z is stronger than that in 
the dot plane allowing one to separate the hole motion 
along z and in the (x,y) plane. As a reasonable sim- 
plification, the potentials acting on each kind of holes, 
n = hh or n = Ih, are modelled by separable functions, 

V n (r) = V z n (z) + Vf(p, lp ), (5) 
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where for convenience we use the cylindrical coordi- 
nate system with the principal axis z and in-plane po- 
lar coordinates p = y^x 2 +y 2 and <p. The eigenstates 
and eigenenergies are found from the time-independent 
Schrodinger equation 

[^r s + V(r)]*(r) = E$(r), 

where ^(r) is a column formed by four envelopes x P m (r) 
with m = -3/2, -1/2, 1/2, 3/2. 

The heavy-light hole mixing is described by the off- 
diagonal components of the Luttinger Hamiltonian Q. 
While calculating the hole g-factor, we take into account 
this mixing in the first-order approximation. In the 
zeroth-order approximation, however, the off-diagonal 
components of T-lr s can be neglected at all, thereby, 
heavy- and light-hole states are completely separated 
in this approximation, so that the four-component hole 
wavefunction ^ has only one non-zero component <J/^ 
with to = —3/2 either 3/2 for n = hh and to = —1/2 
either 1/2 for n = lh^ Due to the assumed separable 
potential d5| the envelopes ^^(r) can be sought in sep- 
arable fornPlMl 



(0) 



Here F[ l (z) describes the hole size-quantization along the 
growth axis, with I = 1,2... is an integer enumerating 
the size-quantized states along this axis, and ipp(p,<p) is 
the in-plane envelope with p = 1,2... labeling the in- 
plane confined states. 

In the further specification of the potential ^ we bear 
in mind that the QDs under consideration have the ap- 
proximate shape of triangular pyramids asymmetric with 
respect to the z — > — z reflection and possess the three- 
fold rotation symmetry in the (x, y) plane, see Fig. [l]and 
Ref. 1181 The asymmetry along the growth axis is mod- 
elled by a triangular potential well, 



-eTz, z > 
-oo, z < 



(7) 



where z = corresponds to the QD base, T is the effec- 
tive electric field breaking the z — > — z mirror symmetry 
and confining a hole along the z-axis, and e = — |e| is the 
electron charge. In this case the z-dependent envelopes 
are expressed via the Airy functions as follows 



Fr{z) = CfM{Zf) , 



(8) 



where C™ is the normalization factor, 



Z? 



Lit 



2m„ iZ eJ" 



1/3 



pi is the l-th root of the equation Ai(— Z) = 0, and the 
effective mass m n ^ z equals to mo/ (71 ± 27) with the ± 
sign corresponding to the light and heavy holes. 

The in-plane potential acting on holes is taken to be an 
analytical function of p and comprise the main parabolic 



terrrPSHSD an( j a cubic correction describing the trigonal 
warping, namely, 

If (^) = -^^(l + /3^cos3^ . (9) 

Here m n 11 = too/ (71 ± 7) are the effective masses for 
heavy (top sign) and light (bottom sign) holes in the QD 
plane, a n are the effective localization radii and the di- 
mensionless parameter j3 characterizes the quantum dot 
triangularity. 

Figure [4] represents equipotential surfaces of the dot 
confining potential defined by Eqs. Q5J) , ([7]), and Q for 
/3 = (a) and /3 = —0.35 (b) and illustrates the transfer 
from the axial to trigonal symmetry. 




HO] (b) 



FIG. 4: Equipotential surface of the dot confining potential 
for the two values of triangularity parameter /3 = (a) and 
/3 = -0.35 (b) 

For the GaAs parameters 71 = 6.98 and 7 = 2.58, one 
has rrihh,z = 0.55to , m lh ^ = 0.08to m hh ^ = 0.1to , 
and TO;/> || = 0.22to . The size-quantization energies for 
the motion along z-axis are given by 



h 2 m 



2m 'n 1 zL 2 l 



1 = 1,2. 



(10) 



It is convenient to introduce the dot effective size in the 
z-direction as L = y/LhhLih- For L — 30 A the values 
of energies of the ground heavy and light-hole states are, 
respectively, w 33 meV and E l J\ sa 63 meV. The en- 
ergies of size quantization in the quantum dot plane are 
labeled as E^ h p and E\h p for heavy and light holes, respec- 
tively. If j3 = the energies of in-plane size-quantization 
form an equidistant set with the spacing 



r, 2 



(ii) 



and the eigenfunctions are those of a two-dimensional 
axially-symmetric harmonic oscillator. In what follows 
we assume that the confining in-plane potential in the 
case of (3 = is identical for the heavy and light holes. 
This implies that TO/^ ya^ = mu l ^af h and the ratio of 
two localization radii for GaAs is ahhl&ih ~ 1.21. Taking 
75 A as a reasonable estimation for ahh we obtain hw^h ~ 
6.5 meV and fvuJih ~ 4.4 meV. These values are small 
as compared to the energy difference E l J\ — E^\ s» 30 
meV of hole confinement in the z-direction, confirming 
the consistency of the model. 
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2. Longitudinal g-factor gni 

The Zeeman interaction of a bulk jT 8 hole with the 
magnetic field B is described by the 4x4 matrix opera- 
tor 51 , resulting for B | [111] in the following expressions 
for the effective g-factorsP^T 



-6x- T?> 



9h2 



2V2q, 



(12a) 



(12b) 



anisotropy of the dot was neglected, i.e., in the frame- 
work of our model we set /3 = and chose the angular en- 
velopes ipp [see Eq. ([6|] as the eigenfunctions of isotropic 
two-dimensional harmonic oscillator. Note, that trigonal 
distortion of the in-plane potential described by a small 
parameter f3 leads to a correction to gni quadratic in 
f3. Indeed, as it follows from the symmetry reasons the 
change of the (3 sign is equivalent to the rotation of the 
triangular QD around z axis by it, in which case ghi re- 
mains unchanged. Estimations show that this quadratic 
correction is negligible. 



where x and q are the dimensionless band structure pa- 
rameters, des cribin g isotropic and cubic contributions to 
Zeeman effect .^Ml In GaAs crystals k = 1.2, q = 0.02P 



It follows from Eq. (12b) that the off-diagonal g-factor 



gya is already present in this model, however, its value 
gh2 ~ 0.056 is too small to describe the experiment. 
More over, similarly to bulk semiconductors and quantum 
wells 1S2153H56] |- ne magnetic field induced mixing of heavy 
and light hole states and quantum confinement strongly 
renormalizes heavy-hole g-factor in quantum dots from 
its value in the bulkp3E31 Indeed, in the presence of the 
magnetic field the wavevector fc in Eq. ^ is replaced by 
k + eA(r)/ch, where A(r) = (1/2)B x r is the vector 
potential of the magnetic fields The vector potential, 
through the off-diagonal components of the matrix Q, 
mixes the ground heavy-hole state with light-hole states. 
Calculation presented in Appendix [A] shows that the lon- 
gitudinal g-factor ghi in the quantum dot is given by 



g hl = -6x - (23/2)g + Ag hl , 



(13) 



where the third term obtained in the second order in L/a 
has the form 



Ag hl = 247 



V 2 +tf(L 2 /2a 2 ) 
, Mi - Mi \/mih,z/m hhiZ 



Here a = a^h, other notations are 



(14) 



d 



m = L I Fr(z)^Fr(z)dz 



& = cf - 



w = / F 1 hh (z)F{ h (z)dz. (15) 



Performing summation in Eq. ( 14 1 over I for the set of 



GaAs parameters 71, 7, mhh ,|| an d toj^ji we obtain for 
9hi 



9 hi 



-6x- (23/2)^ + 4.5 + 6.7 



L 2 



(16) 



where the third and fourth coefficients are calculated 
numerically. For the g^i calculations the in-plane 



3. Off-diagonal longitudinal g-factor g^i 

The pyramidal shape of the quantum dot with the 
three-fold rotation axis z results in the magnetic-field in- 
duced mixing of heavy holes with opposite angular mo- 
mentum projections ±3/2. The trigonal distortion of 
the in-plane parabolic confinement described by the term 
oc f3p 3 cos 3ip in Eq. |9| leads to the lowering of symmetry 
of the ground heavy-hole state envelope function, which 
can be represented in the lowest order in |/3| <C 1 as 



?Pi h oc exp 



2a 2 ) 



l+P^C(p) cos3</> 
a 6 



(17) 



where function C(p) is finite at p — and only terms 
up to linear order in (3 are retained. The variational 
procedure with the space-independent C{p) = Co gives 
for Co a value of —1/6. 

By means of the second-order perturbation theory we 
obtain for the magneto-induced mixing of 3/2 and —3/2 
heavy-hole states [cf. Eq. Q]: 



^B,3/2,-3/2 = 2^ B 9h2B z 



(18) 



E 

lp± 



1 ±l/2;Z 



I -3/2;ll 



rphh 



E: 



11, 



We note that while replacing in rir s the hole wavevector 
k by k + eA/hc only linear in A terms should be retained. 

The details of evaluation of the off-diagonal parameter 
gh2 are presented in Appendix [5] Similarly to the cal- 
culation of g^i, in the limit L/a ^ 1 one can neglect 
the dependence of hole energy on the in-plane quantum 
numbers p with the result 



9h2 



18^7 



2^,2 



V1C1 



7TLqL 



I 



E: 



in 



rphh 



(x - iy) 



d_ 



dy 



(19) 



This equation clearly shows that the value of gy^ is 
nonzero provided only that (i) the dot base has a tri- 
angular shape, i.e., (3^0, otherwise integration over the 
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in-plane coordinates vanishes, and (ii) the dot is asym- 
metric in the z direction, otherwise the product rjiQ van- 
ishes due to parity. 



Using the variational function (17) with C(p) = —1/6 
we obtain 



9h2 = 36/37 



m hh, z mf hz L 



m 



E 



m - Ml ^/mih, z /m hh ^ 



(20) 



The numerical summation in Eq. ( 20 1 yields 
g h2 « 10.2/3- , 



(21) 



where the small contribution 2\/2q is neglected. For 
reasonable values j3 = —0.2, L/a — 0.3 we obtain 
\ghi \ ~ 0.62, in good agreement with experiment. 



IV. DISCUSSION 



Figure [5] displays the g- factors ghi and gh 2 calculated 
as a function of the ratio L/a according to Eqs. ( 16 1 and 
(21 1. The figure shows that the values of \gh2\ vary in the 



range of 0.3-^0.7 for a reasonable range of the parameter 
/3 and the ratio L/a in good agreement with the experi- 
mental observations. The diagonal g- factor, g^i, changes 
its sign at some critical value of L/a which depends on 
the material parameters and quantum dot shape. We 
note that for the particular form of the Zeeman Hamilto- 
man, Eq. p}, the sign of ghi cannot directly be extracted 
from the PL data in the magnetic field. Relying on our 
estimations we believe, however, that the value of ghi for 
studied dots with small values of L/a is negative. 

Accurate numerical estimations of effective g-factors 
require the knowledge of spacing between size-quantized 
levels of heavy and light holes in the z-direction. 
The value of this spacing can be estimated using the 
anisotropy of the electron g-factor, i.e., the difference be- 
tween its longitudinal (g ej ||) and transverse (g e .±) com- 
ponents. It is well known that the g- factor of an e lectron 
confined in a quantum dot is strongly anisotropic] 44 * 65 * 6 -^ 
and the value of that anisotropy is directly related to 
the energy spacing between heavy- and light-hole states. 
Taking for the experimental values of g e< \\ and g e ^± their 
relative difference (g e ,i_ — g e ,\\)/9e,\\ ~ 40% and using for- 
mulae of Ref. UH the energy spacing E l z \ — E z \ can be 
estimated as 20 30 meV, which is in a good agreement 
with the values obtained in our simple model of the con- 
finement in the z-direction with L ranging from 25 to 
35 A. 



It should be noted that Eqs. (16 1 and (211 are derived 
assuming L/a <C 1. We presume that they can reason- 
ably be extrapolated up to QDs with L/a of the order of 



unity and use the following working equations 



I? 



g h i(L/a) = g° hl + c 1 - 
a z 



g h2 {L/a) = c 2 j3 



L 



(22) 



where gl x — <?/ii(0) and the coefficients c\,c 2 are con- 
sidered as fitting parameters. The absence of zero-order 
term in the dependence of gh 2 on L/a is expected since, 
in the band model under consideration, the heavy-hole 
mixing in quantum wells is negligible. In contrast, the 
diagonal component ghi is strongly renormalized even 
in quantum wells. The quadratic term in (16 1 appears 

Eq. ((All 



due to the small term _En 1 '? — E l n h in the denominators of 
_ _ _ _ II ,p 

Eliminating the variable L/a we directly relate 
the diagonal and off-diagonal components of the ^-factor 
as 



9h2 = c 2 (3\ 



! 9hi 



Cl 



(23) 



The product c 2 f3/^/c\ can be considered an indivisible 
parameter. We, however, take c\ = 6.7, c 2 — 10.2 as 
follows from Eqs. ( fl6| , (21 ) and fit the value of (3 to the 
experimental data. Equation (23 1 predicts a particular 



correlation between the diagonal and off-diagonal hole g- 
factors. Figure [6] represents experimental results as well 
as theoretical fit with parameters introduced in the cap- 
tion. As mentioned above, the experimental values of 
ghi are taken with negative signs. One can see that the 
theory well reproduces the experimental relation between 
g h2 and g h i- 




P = -0.2 



-0.15 



-0.1 



FIG. 5: The off-diagonal <;-factor gh2 as a function of L/a 
depicted for three values of the trigonal parameter f3: —0.1, 
—0.15 and —0.2. The inset shows the behavior of ghi- 

The systematic differences between the off-diagonal g- 
f actors gh 2 for X + and X~ trions seen in Table |l] indicate 
a non-negligible role of the Coulomb interaction. Indeed, 
the value of ghi found from the spectrum of X + trion 
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FIG. 6: Correlation between the values of g^i and \gh2\- Sym- 
bols represent experimental data for different GaAs dots mea- 
sured on the neutral as well as on charged complexes, dashed 
curve stands for theoretical fit with parameters j3 — —0.15, 
= —1.5, ci = 6.7, C2 = 10.2. Thin lines are calculated 
for the same set of parameters, except for /3 = —0.13 and 
—0.17. The corresponding values of L/a (for /3 = —0.15) are 
presented at the right-hand axis. 



radiative recombination corresponds to the off-diagonal 
g-factor of a resident hole, while the value gya extracted 
from the X~ spectrum gives the g-factor of the hole in 
a trion. However, the detailed investigation of Coulomb- 
driven effects is beyond the scope of the present paper. 

As far as we use Hamiltonian Eq. Q in the spherical 
approximation, the [111] growth direction is nominally 
equivalent to [001] direction, except for natural trigonal 
shape of the c onfining potential of QDs grown in the for- 
mer cas e! 18 * 21 !. A nonzero value of ghi is also obtained if 
we neglect the QD triangularity but make allowance for a 
cubic anisotropy of the Luttinger Hamiltonian, described 
by the nonvanishing difference 5j — 72 — 73 which is a 
small parameter for majority of semiconductor materials. 
Transformation of this Hamiltonian from the crystallo- 
graphic axes [100], [010], [001] to the frame x,y,z results 
in its matrix form M| where H and / acquire corrections 



Therefore, we can conclude that a sizable off-diagonal 
heavy-hole g factor, gh2, indeed results from the quan- 
tum dot shape in the form of triangular pyramid. Note 
that this is independently confirmed by an absence of 
the heavy-hole mixing in [111] grown disk-shaped dots 
studied experimentally in Ref. 1671 



V. CONCLUSIONS 

In this paper we have shown that GaAs droplet dots 
grown along [111] axis demonstrate striking magneto- 
photoluminescence spectra. In longitudinal field B \ 
[111], each electron-hole complex, a neutral exciton X° 
and positively X + or negatively X~ charged exciton, 
shows four emission lines with definite circular polariza- 
tions, instead of two lines for commonly studied (001) 
dots. Moreover, in transverse field B _L [111] only two 
linearly polarized lines are observed. These effects are 
described by the magneto-induced mixing of the heavy 
holes for B | [111] and the absence of such mixing for 
the transverse field. 

A microscopic theory of the mixing has been devel- 
oped in the framework of Luttinger Hamiltonian ap- 
proach adopted to describe the hole quantization in the 
dots with a shape of triangular pyramid. Our theory is 
in good quantitative agreement with experimental data. 
We have established that it is the shape of the dot which 
results in off-diagonal components of the g-factor tensor 
and gives rise to the novel magneto-optical phenomena 
in droplet (111) quantum dots. 
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SH 



SI = 



h 2 Sj 

h 2 6>y 



(k 2 _ - 2V2k z k+^j , 



(24) 



where k± — k x ± ik y . These corrections lead to nonzero 
values of gya even for the disk- and lens-shaped dots as 
well as for the [111] grown quantum wells. However, the 
corresponding contribution Sgya oc Sj is indeed small 
and, in the case of trigonal dots, constitutes a fraction 
< 10% of that obtained from Eq. ( 20 I . As mentioned 
above, an additional mechanism of the magneto-induced 
mixing of the ±3/2 hole states suggested in Ref. I2~4l and 
based on a cubic-fc contribution to the hole Hamiltonian 
also gives too small value of 1 ■ 



Appendix A: Calculation of gni 



The additional term Ag^i in Eq. (13 1 for the diagonal 



g-factor g^i is found in the second order of perturbation 
theory from the following relation 



1 



-Ag hl (i B B z 



(Al) 



E 

i P ± 



*3/2;lll^I 



I ^±l/2;Zp 



ylh 

*±l/2;Zp 



I 3/2;ll 



E 



iii 



E 



lh 



where the summation includes all the light hole states 
described by the two indices I and p. Using the hole 
envelopes one can rewrite (Al I in the form 
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Ag hl = {2/fi B B z 



i P 



hh 
1 



E 



lh 



phh 



■n 



H\Fl h 



lh 



F IK 



i>l h \H' 



E 



lh 



hh 



F? n i/j? n \I\Fl h 4> 



(A2) 



Ff h ip lh \I*\F? n ^1 



hh 



Here H and / are the components of the Hamiltonian (|4| 
modified in the presence of magnetic field: 



H = -V3Bk z 



ky H fel-^X _ lAy) 



ch 



2e 



ik y ) 2 + ^ {(k x - ik y ){A x - iA y )} 



k a should be replaced by the differential operators 
—id/dx a , the curly brackets denote a symmetrized prod- 
uct of the corresponding operators. In the calculation of 
the sum (A2| one has to retain only linear contributions 
in A. In the limit L/a « 1 the energy denominator in 
Eq. (A2| can be expanded in powers of this ratio. Per- 



forming summation over the index p we directly arrive at 
Eq. (14) from the main text. 



Appendix B: Calculation of gu2 



arrive at 



9h2 



m L ^ EV\ + Eff h , - E lh , - E" 1 
' lp z . 1 ||,i 2 ,' II, v 



i v 

4>i h Iv 2 . 



J,hh\ 



+ (^|V_K)(<|{V_r_}|^)) , (Bl) 



where r_ = x — xy oc e lip , V_ 
V 2 _ (x (T 2 ^. Note that at /3 



= d/dx — id/dy oc e~ lv , 
the mixing described 
by Eq. (Bl ) vanishes because in this case the matrix ele- 
ments of operators proportional to e~ ltp and e~ 2lip cannot 
be nonzero simultaneously. 

Retaining the first-order term in the expansion of 
in powers of L/a -C 1 we can neglect t he in-plane en- 
ergies Eu in the denominators of Eq. (Bl). Then the 



Starting from Eq. ([18]) and using the components of summation in Eq. (Bl I over the index p canbe performed 



the Luttinger Hamiltonian H,I in the magnetic field we analytically, and we directly arrive at Eq. (19) 
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